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Abstract 

Vortex structure of d^2_y2-wave superconductors is microscopically analyzed 
in the framework of the quasi-classical Eilenberger equations. If the pairing 
interaction contains an s-wave (d^j^-wave) component in addition to a (1^2 _y2- 
wave component, the s-wave {dxy-wave) component of the order parameter 
is necessarily induced around a vortex in (i^.2_j^2-wave superconductors. The 
spatial distribution of the induced s-wave and dxy- wave components is cal- 
culated. The s-wave component has opposite winding number around vortex 
near the d^2_j^2 -vortex core and its amplitude has the shape of a four-lobe 
clover. The amplitude of daj^-component has the shape of an octofoil. These 
are consistent with results based on the GL theory. 
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In the last years, a number of investigations were carried out theoretically and experimen- 
tally to identify the symmetry of pairing state in high-T,, superconductors. Although precise 
pairing symmetry has not been determined yet, it is recognized that c?a,2_y2-wave symmetry 
is most probable.0 Recently, the vortex structure of the rf-wave superconductors attracts 
much attention because it may have different structure from that of conventional s-wave su- 
perconductors. One of the related topics is a possibility that other components of the order 
parameter may be induced around a vortex in (i-wave superconductors. Based on symme- 
try considerations, Volovik suggested this possibility.! In (ia;2_j,2-wave superconductors, the 
amplitude of s-wave pairing should be contained at the core region of the vortex, since the 
(ia;2_y2-wave vortex has the same symmetry as that of the opposite winding s-wave compo- 
nent. Solving the tight-binding Bogoliubov-de Gennes equation self-consistently, Soininen, 
Kallin and Berlinsky calculated the vortex structure on a 16x16 lattice and showed that 
s-wave component with an opposite winding number is contained around the vortex core 
in (i2.2_j^2-wave superconductors.il Ren, Xu and Ting derived the two components Ginzburg 
Landau (GL) equations for s- and (i-wave superconductivity from the Gor'kov equations, 
and studied the mixing of the s-wave component near a vortex in d-wave superconductors.^ 
They suggested that the winding number of the induced s-wave component around the vor- 
tex is 3 far from the vortex core and —1 near the center of vortex, and its profile has the 
shape of a four-lobe clover. 

According to the general consideration based on the GL theory, it is possible that the s- 
wave component is coupled with the d-wave component through the gradient terms. There- 
fore, s-wave component is induced wherever the (i-wave order parameter spatially varies, 
such as near the vortex. The mixing due to the same mechanism also occurs near a surface 
of (i-wave superconductors, where ci-wave component decreases on approaching the surface 
and s-wave component is induced in the surface region. This surface effect was considered 
by Matsumoto and Shiba using the self-consistent quasi-classical Green function formalism.i 
They also considered the mixing of d^y-wave component in addition to the mixing of s-wave 
component. 
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The purpose of this paper is to analyze s- and d^y-wave components induced around 
a vortex in (ij.2_y2-wave superconductors using the quasi-classical Eilenberger equationsjl 
which can be applied at arbitrary temperatures. The quasi-classical calculations on the 
vortex structure were carried out for conventional s-wave superconductors by Pesch and 
Kramer ,0 and Klein,! and for pure dx2_y2-wave superconductors by Schopohl and Maki,ill2l 
and the current authors.0 Here we consider the case of an isolated vortex under a mag- 
netic field applied parallel to the c-axis (or ^-axis) in the clean limit. The Fermi surface 
is assumed to be two-dimensional, which is appropriate to high- Tc superconductors, and 
isotropic for simplicity. Throughout the paper, lengths and energies are measured in units 
of the coherence length ^ and the uniform gap Aq at T = 0, respectively. 

First, we solve the quasi-classical Eilenberger equations to obtain the Green functions 
for a given pair potential.0 For the dx'^-y^ symmetry and not too low temperatures, the pair 
potential may be assumed by 

A(^^,r) = A(0,r)e*^ 

A(^, r) = A(T) tanh(r) cos(2^), (1) 



where r = \/x^ ^ is the distance from the center of vortex line, Q is the angle of k-vector 
with the a-axis (or x-axis), and the phase of the pair potential around vortex center is 
given by e**^ = {x + iy)/r. The quasi-classical Green functions with the Matsubara frequency 
ujn = (2n+l)7rT are obtained by solving the Eilenberger equations, which are given as follows 
in the gauge where pair potential is real, 

{u;n + ^{d\\+zd\\cf>)}f{u;n,9,r) = A{9,r)g{u;n,e,r), (2) 
- Ud\\ - zd\\(p)}p{u;r,,9,r) = A{9,r)g{u;„.,9,r), (3) 



g{u;n,9,r) = (l- f{u;^,9,r)p{u;^,9,r)y^\ (4) 

where Iieg{uJn,9,r) > and d\\(j) = —r±/r'^. Here, we have taken the coordinate system: 
u = cos 6'x + sin ^y, v = — sin ^ic + cos 6'y, thus a point r = xx + yy is denoted as r = 



r|jU + r_LV. The anomalous Green functions / and /''' in Eqs.(||)-(^ are related to the usual 
notations / and as / = /e*"^ and = /^e"*"^. 

Instead of solving Eqs.(0)-(^), it is more convenient to use the following parameterization 
devised by Schopohl and Maki,i 

1 + aF ^ 1 + ab' 



From Eq.(^), g is given by 



1 - ab 



Substituting Eqs.® and (P) into Eqs.(||) and (H), we obtain the Riccati equation for a: 

di\a{ujn, e, r) - A(^, r) + {2cu„ + id\\(j) + A{e, r)aK, 6, r)}d{un, 6, r) = 0. (7) 

The other unknown quantity b is related to a by symmetry: b{r\\) = a{—r\\). Far from the 
vortex core, a is reduced to the value of a spatially homogeneous situation without magnetic 
field: 



aoo{e)= ' " L' 1 {Un>0). (8) 

To obtain the quasi-classical Green functions, we integrate Eq.(|^) along the trajectory where 
r_L is held constant, using Eq.(H) as the initial value. 

Second, we calculate the pair potential from the resulting quasi-classical Green functions 
by the self-consistent condition 

/-27r f]f) 

A{e, r) = No2nT ^ / ^V{e, 0')/^, 0\ r), (9) 

where Nq is the density of states at the Fermi surface. The pair potential and the pairing 
interaction are decomposed into s-, dx2_y2- and dxy-vf&ve components, 

V{e, 9') =Vs + V^2_y2 cos(2^) cos(2^') + V^y sm{29) sin(2^')> (10) 
A(e, r) = A,(r) + A,2_y2{v) cos{2d) + A,.,(r) sin(2e). (11) 
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Substituting Eqs.(|TO|) and (|TTD into Eq.@, we obtain 

de 

Wn>0 



A,(r) = KiVo27rT ^ / — />„, 6, r), (12) 



/■2tt _ 

A,2^y2{r) = V,2_y2No2nT —f{un,e,r)cos{29), (13) 



UJn>0 



Ky{v) = V,yNo2nT / — />„, 6, r) sin(2^^). (14) 



UJn>0 



As already mentioned, at not too low temperatures, T/Tc > 0.5, the resulting profile of 
A^.2_j^2(r) calculated from Eq. (|1^) after solving the Eilenberger equations is almost the 
same as that of Eq. (0), ensuring the self-consistency in that temperature region. 

As seen from Eqs.(|T2|) and (0), the Green function / solved under the given pair poten- 
tial immediately yields the induced s-wave (or dxy-wave) component. Figures |I| and ^ show 
the s-wave component As(r) induced around a vortex in dx2^y2-wave superconductors. As 
seen from Fig. |l|, the amplitude has the shape of a four-lobe clover. Near the vortex center 
\As\ oc r and far from the vortex core \As\ oc r~^. As shown in Fig. ^, the term with e"^**^ 
is dominant near the vortex center and the term with e^**^ is dominant far from the vertex 
core. The induced s-wave component, therefore, can be written as 

A,(r) = A,(r)e*'^ = (^Ci{r)e~^'^ + C2{r)e^'''')e'^, (15) 

where ci(r) and C2(r) are factors depending on r. Near the center of vortex |ci| > |c2|, and 
far from the vortex core |ci| < \c2\- These are consistent with the results given by Ren, 
Xu and Ting based on the GL theory.Blli On lowering temperature, the amplitude |A5(r)| 
increases and the inner area where e"^**^ is dominant spreads out. 

Figures | and || show the d^y-wave component A^y(r) induced around a vortex in dx2_y2- 
wave superconductors. As shown in Fig. ^ the amplitude has the shape of an octofoil. Near 
the vortex center |As| oc and far from the vortex core |As| oc r^^. As seen from Fig. ^, 
the term with e"^**^ is dominant near the vortex center and the term with e^*"^ is dominant 
far from the vortex core. The induced dxy-"wave component, therefore, can be written as 



A,,(r) = Ky{v)e"t' = (d,{r)e-'"t' + d2{r)e^'^)e't (16) 

where di{r) and d2{r) are factors depending on r. Near the center of vortex \di\ > \d2\, 
and far from the vortex core \di\ < \d2\. Our resuhs of the induced (ia-j^-component are also 
explained by the GL theory, if the non-local correction terms are included. As far as the 
gradient terms are concerned, the mixing of d^2__y2- and d^^-components is absent in the 
usual second derivative terms and first appears in the fourth-order derivative terms. 

Let us now interpret our microscopic calculations so far in terms of the GL framework for 
ease of the understanding. The GL equations are derived from the Gor'kov equations by the 
same calculation as in the case of s- and (ij:2_j^2-components.S For d^^^yi- and (i^^jy-components 
which are simultaneously non-vanishing, we obtain 

+ |A,.J2A,, + ^|A,.2_,2|2a,^ + ^^l2_y2^ly = 0, (17) 

- In (^) A.._,. - {dl + 9J)A..._,. - 7{ {^-{dl + dlf - dldi)^.^.^y. + \d^dy{dl - dl)Ky 

2 1 

+ |A^.2_j^2| A^.2_y2 + -\^xy\ A2.2_y2 + -A^,j^A*2_y2 = 

in the dimensionless form, where 7 = 62C(5)/49C(3)^, axy = 2{(yxyNo)~^ — {Vxi-y^No)'^}. 
Far from the vortex core, we can assume A^.2_j^2 = AqoC*"^, and the leading order terms of 
the GL equations (0) and (|18D are written as 

(axy + ^)Axy + ^e'^'^Al^ - ^d.dyidl - d'^^ = 0, (19) 

-ln(T,/T)Aoo + AL = 0. (20) 

Solving Eqs. ([ToD and (^^, we obtain 

Axyiv) = -^(rfie-^^^ + 4e^^^)e^^ (21) 

where d[ = d'ia.y + 3A^), 4 = rf'(7a., + SA^), d' = {{a^y + |A^)2 - {lAm-'^A^ 
and Aqo = {ln(Tc/T)}-'^/^. Since jcigl > \d[\, the term with e^**^ is dominant. Near the center 



of vortex, the leading order of Eqs. ([T7|) and ([Tq) are second and fourth order derivative 
terms, and the general solution of the GL equations is given as follows. For a d^2_y2-wme 
component, l^x'^-y^ = {(^o^ + ^o'^^ + 0(r^))e*'^, where qq and 60 are constants. For a dxy- 
wave component, assuming the form of Eq. ( |T^ ) for Axy, we obtain di = d^r^ + O(r^) 
and d2 = O(r^), where d^ is a constant. Therefore, the term with e"'^*'^ is dominant and 
\Axy \ oc r^. These results are consistent with our quasi-classical calculations. 

Our numerical results are quantitatively valid for small \ Vs/Vx2-y2 \ and |K.j//K:2_y2 1, and 
for not too low temperatures. From Eqs. (|12|) and ([l^), As(r) and A^,j^(r) are proportional 
to Vs and Vx^;, respectively. If |Kj/V^2_j^2 1, |\4y/K'2-2;2| ^ 1, ^^(r) and A^y^r) are negligibly 
small. The s- and d^y-wave components, then, do not affect (i^2_y2-wave superconductivity. 
In the case that Vs or Vxy is comparable to V^2_y2, self-consistent calculations including s- 
and dxy-wave components are needed. It should be noticed that at low temperatures, a 
self-consistent calculation for Ax'2-y2[r) is needed since the pair potential may deviate from 
Eq. (g). 

A vortex in dxi-yi-v/ave superconductors has a four-folded symmetric structure by the 
mixing of the s- or (i^j^-wave component. However, this anisotropy clearly appears only 
when Vs or Vxy is comparable to Vx2^y2. It may not be probable that this condition is 
satisfied in high-Tc superconductors. Lastly, we note that the cylindrical symmetry around 
a vortex line is spontaneously broken in (ij.2_j^2-wave superconductors even when Vs and Vxy 
are absent. The four-folded symmetric structure of (i^.2_^2-wave vortex becomes clear on 
lowering temperature, which is confirmed by our quasi-classical calculations.0 

We would like to thank Professor T. Ohmi for valuable discussions. 
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FIGURES 

FIG. 1. (a) Amplitude of the s-wave component, |As(r)|/(T4/T42_^2), induced around a vortex 
in dx2-y2-wave superconductors at T/Tc = 0.5. It has the shape of a four-lobe clover, (b) The core 
region of (a) is focused. Near the center of vortex, |As| oar. 

FIG. 2. Phase of the induced s-wave component, argA5(r). Far from the vortex core, the 
term with e^'*^ is dominant. Near the center of vortex, the term with e"^**^ is dominant. 

FIG. 3. (a) Amplitude of the dxy-wawe component, \Axy{r)\/{Vxy/Vx2_y2), induced around a 
vortex in dx2_y2-wave superconductors at T/Tc = 0.5. It has the shape of an octofoil. (b) The core 
region of (a) is focused. Near the center of vortex, \Axy\ oc r^. 

FIG. 4. Phase of the induced dxy-w&ye component, aTgAxy{r). Far from the vortex core, the 
term with e^^'^ is dominant. Near the center of vortex, the term with e"'^^'^ is dominant. 
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Fig. 3 (a) 



Fig. 3 (b) 



Fig. 4 



